Turbulent Relative Dispersion in Two-Dimensional Free Convection Turbulence 
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The relative dispersion process in two-dimensional free convection turbulence is investigated by 
direct numerical simulation. In the inertial range, the growth of relative separation, r, is expected 
as {r^{t)) oc t^ according to the Bolgiano-Obukhov scaling. The result supporting the scaling is 
obtained with exit-time statistics. Detailed investigation of exit-time PDF shows that the PDF is 
divided into two regions, the Region-I and -II, reflecting two types of separating processes: persistent 
expansion and random transitions between expansion and compression of relative separation. This 
is consistent with the physical picture of the self-similar telegraph model. In addition, a method 
for estimating the parameters of the model are presented. Comparing two turbulence cases, two- 
dimensional free convection and inverse cascade turbulence, the relation between the drift term of 
the model and nature of coherent structures is discussed. 



I. INTRODUCTION 

Relative dispersion of passive particles in turbulent 
flows is one of the fundamental problems in turbulence 
research. It characterizes the transport and mixing prop- 
erties of turbulence and is important from both theoreti- 
cal and practical points of view. Reflecting universal be- 
havior of turbulent fluctuations, relative dispersion also 
has some universalproperties because of its locality-in- 
scale nature In particular, the dispersion pro- 

cess exhibits anomalous dispersion in the inertial range. 
This is first observed by Richardson (1926) 0, and since 
then, a number of theoretical, experimental, and numer- 
ical investigations have been devoted to understand and 
model relative dispersion process jsl. ^G*. "t*. "q UioIIiiIIT^ . 
However, comprehensive understanding has not been ob- 
tained yet. 

Recently, a few works focusing on underlying mecha- 
nism of the anomalous dispersion were reported. In the 
inertial range, the mean free path, l{r), the mean length 
for persistent expansion of relative separation without 
changing its moving direction, is an order of relative sep- 
aration itself, r: l{r) = Pgr ^^uM- Here, Pg is a non- 
dimensional constant called the persistent parameter. In 
two-dimensional inverse cascade turbulence (2D-IC), Pg 
is estimated as 0.87 This means separating motions 
are not purely diffusive but composed of an appreciable 
amount of persistent motions. In addition, it was also re- 
ported that there is a relation between stagnation-point 
structures and Richardson's law and that disper- 
sion process is described by persistent streamline topol- 
ogy *lu\. From these results, it is expected that coher- 
ence of turbulent field, which must share its origin with 
fine coherent vortical structures such as worms in three- 
dimensional Navier-Stokes (3D-NS) turbulence, has a sig- 
nificant role in turbulent relative dispersion. 

The correlations in turbulence are characterized in 
scale-space due to their self-similarity, and are not made 
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disappear by coarse graining in real space. Reflecting 
these correlations, relative separation moves persistently 
to some extent, so that, unlike the Brownian motion, rel- 
ative separation process should not be described only by 
random collision motions even as an approximation 13J. 
In other words, the characteristic length cannot be de- 
fined because the mean free path, l{r), varies depending 
on the spatial scale. Whereas there are several exper- 
iments and numerical simulations of which results are 
rather close to the prediction of Richardson's diffusion 
equation 0,11,^3 ^^at closely relates to random coUision 
motions. Thus, these results raise a question; how are 
the effects of persistent motions wiped out? In the pre- 
vious paper |l6j . we have introduced a self-similar tele- 
graph model of turbulent relative dispersion, and showed 
that the separation PDF can be close to the prediction 
of Richardson's diffusion equation for slowly separating 
particle pairs even in the presence of persistent motions. 

In the present paper, we check the consistency of the 
physical picture of the self-similar telegraph model by 
carrying out direct numerical simulations (DNS) of 2D 
free convection (2D-FC) turbulence instead of 3D-NS tur- 
bulence. This is because (i) DNS of 3D-NS turbulence 
requires extremely large computer resources, so that it 
is difficult to track particles for a long time, which is 
necessary to investigate dynamical properties of rela- 
tive dispersion process, and (ii) 2D-FC turbulence has 
both statistical and dynamical characteristics similar to 
those of 3D-NS turbulence (see Appendix A) 
Among them, the existence of coherent structures, which 
are approximated by the Burgers T- Vortex layer, is no- 
table. This is a crucial difference from coherent struc- 
tures in 2D-IC turbulence that are nested cat's eye vor- 
tices Pldj. Therefore, comparing the results of the 2D- 
FC case with those of the 2D-IC case, we can investigate 
the effects of coherent structure on turbulent relative dis- 
persion. This comparison gives a physical meaning of the 
drift term of the self-similar telegraph model. 

The inertial range achieved by our DNS is not so wide 
that the relative motions of particle-pairs in the dissipa- 
tion, the inertial, and the energy containing scales are 
not sufficiently resolved by usual fixed time statistics. In 
order to investigate the scaling natures of relative dis- 
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persion in such a narrow and limited inertial range, we 
utilize exit-time statistics introduced into research of tur- 
bulent relative dispersion by Boffetta and Sokolov (see 
Appendix B) "s", Ty, ^2^. By detailed investigation of the 
PDF of exit-time, wc show that the PDF is divided into 
two region, the Region-I and -II, corresponding to persis- 
tent expansion and random transition between expansion 
and compression of relative separation, respectively. This 
result agrees with the picture of the self-similar telegraph 
model. In addition, we provide a method for estimation 
of the parameters of the self-similar telegraph model by 
using exit-time PDF. 

In the following sections, first we provide a brief re- 
view of the self-similar telegraph model in §2. Section 
3 presents a summary of the numerical scheme and pa- 
rameters of our DNS of 2D-FC turbulence. Some of the 
results by the fixed-time and the exit-time statistics are 
provided and discussed in §4. Concluding remarks are 
made in §5. In addition, some properties of 2D-FC turbu- 
lence and exit-time statistics are presented in Appendix 
A and B, respectively. 



II. SELF-SIMILAR TELEGRAPH MODEL AND 
PALM'S EQUATION 

In the previous paper |16|. we have introduced a self- 
similar telegraph model for turbulent relative dispersion, 
which is a model describing the evolution of the PDF of 
relative separation of particle pairs in the inertial range. 
The relative separation of two particles, r{t) is defined as 
follows: 



r{t) = \xi{t) - X2{t)\, 



(1) 



where Xi{t) and X2{t) are the Lagrangian positions of 
the particles. The model is based on Sokolov's model 
and consists of persistent expansion and compression of 
relative separation, r, according to the relative velocity, 
v{r) = Ar^~'', where A and s are a dimensional con- 
stant and a scaling exponent, respectively: s — 2/3 for 
Kolmogorov scaling and s = 2/5 for Bolgiano-Obukhov 
scaling. The transition rate from expansion to com- 
pression and that from compression to expansion are 
given by X'^ /Tc{r) and X~/Tc{r), respectively. Here 
Tc{r) is a characteristic time scale at a spatial scale r, 
Tc{r) — A^^r", and are the inverses of persistent 
parameters, P_|_ for expansion and P_ for compression. 
Introducing A = A+ + A~ and (5 = A+ — A^, the evolu- 
tion equation of the separation PDF, P{r,t), is derived 
as follows: 



A df^ 
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-o— \v[r)P\ , (2) 
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where T^ir) is 
A\~^r'^~\ and a 



Richardson's diffusion coefficient, 
= (d - 2s -f = ctA-1. The 



parameters of the model are A and 8: A represents the 
strength of persistency of moving directions, and b does 
the difference in persistency between expansion and 
compression. Therefore, A and (5 reflect the strength and 
structure of coherence of the flow, respectively. 

For slowly-separating particle pairs, i.e., in the case of 
r <C (r^)^/^, the flrst term of the l.h.s. of Eq. |(2Jl can be 
neglected and the approximated equation is given by 



dP 
'dt 



d_ 

dr 



D{r) 



d_ 

dr 



P 



d 

r 

dr 



[v{r)P]. (3) 



This form of the equation was first derived by Palm 
and is also the same as the diffusion equation of Goto- 
Vassilicos model 0. We call Eq. © Palm's equation. 
The similarity solution of Eq. © is 



P{r,t)^Cpt- 



Xr" 



exp 



1 Xr' 
At 



(4) 



where Cp is the normalization factor. 

Because the tail of the exit-time PDF, PE{TE;r, p), 
consists of slowly-separating particle pairs, it is calcu- 



lated from Palm's equation 
by Boffetta and Sokolov 

PE{TE;r, p) is given by 



lO with the method used 
. The asymptotic form of 



PsiTEir, p) 
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Te 
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(5) 



where ji^^„ is the n-th zero of the i/-th order Bessel func- 
tion and {TE{r; p)) is the mean exit-time from r to pr: 



{TE{r-p)) 



X 



1 



(p«-l)r^ 



(6) 



As{2s-5) 

Note that the mean exit-time calculated from the self- 
similar telegraph model, Eq. is the same as Eq. 10. 
This is because the mean exit-time is calculated from a 
steady solution of the equation 8], and the solution of 
Eq. 10) is the same as that of Eq. By comparing the 
tail of the exit-time PDF obtained by DNS and Eq. |(SJ), 
we can estimate the value of 5. 

The last term of the r.h.s. of Eq. ||2J) is a drift term 
consistent with the scaling law, and the direction of the 
drift is determined by a. The parameter a consists of 
two parts, the "scaling-determined" one, d— 2s, and the 
"dynamics-determined" one, S. In order for Eq. (j2Jl to 
recover Richardson's equation, the drift term has to dis- 
appear, which means S = 2s — d = Sq. We call this 
case the Richardson case or the zero-drift case, where 
the parameters of the model reduce to one, A+; A~ is 
determined by the relation A^ = X^ — Sq = Xq . 



III. NUMERICAL SIMULATION 

In this section, we explain the method of DNS used in 
the present work and show basic properties of turbulent 
fields produced by our simulation. 
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FIG. 1: Rescaled energy and entropy spectrum E{k) and 
S{k) obtained by DNS at resolution N = 1024 (NVIO) 
and iV = 2048 (NVU). The straight lines refer to the 
Bolgiano-Obukhov scaling: E(krie) = rig^r^'^ E{k), S{ki]e) = 
Ve'^e'r,-'S{k). 



TABLE I; Parameters used in the present DNS. 



label 


TV 


u 


/o 


do 


kd 


NVIO 


1024 


1 X 10"" 


5 X 10"^ 


1 X 10"^ 


3 


NVll 


2048 


5 X 10"^ 


5 X 10"^ 


1.5 X 10"^ 


10 



We generate turbulent field by DNS of the vorticity 
equation with a large-scale friction term F^, and the tem- 
perature equation with a large-scale forcing term Ff. 



dT 



vAlu - 



dT ^ 

+ Fd, 

ox 



— + iu-\7)T = nAT + Ff, 



(7) 

(8) 

(9) 



where w, T, u, v, k, a, and g represent the vorticity 
field, the temperature field, the velocity field, the kine- 
matic viscosity, the thermal diffusivity, the thermal ex- 
pansion coefficient, and the gravitational acceleration, re- 
spectively. The large-scale forcing term used here is 



Ff{x) = 4/oCos(2a;)cos(2j/), 



(10) 



where /o is a constant. The large-scale friction term is 
written in the Fourier space as 



Fd{k) 



'-^u;{k) {0<\k\<kd), 
(otherwise). 



(11) 



where do, k, Fd{k), and u{k) are a constant, the wave 
number vector, the Fourier mode of the friction term, 
and the Fourier mode of the vorticity field, respectively. 
Our DNS is performed on a 27r x 27r domain with the 



doubly periodic boundary conditions at resolutions iV^: 
N = 1024 (NVIO) and 2048 (NVll). We employ a 
pseudo-spectral method for accurately calculating con- 
volutions and spatial derivatives, and 4-th order Runge- 
Kutta method for time integration. Aliasing error is re- 
moved by adopting the phase-shift method (NVIO) and 
the 3/2 method (NVll). All results presented in this pa- 
per are obtained for statistically stationary and locally 
isotropic turbulence. We summarize the parameters used 
in our simulation in Table and characteristic quantities 
of generated turbulence in Table ITU 

Figure ^ shows the entropy and energy spectra ob- 
tained by our DNS. The spectra are rescaled with entropy 
dissipation scales and the entropy dissipation rate, and 
there is a region consistent with the Bolgiano-Obukhov 
scaling, Eqs. (|A5|) and HA6|I . (see Appendix A). We call 
this region the inertial range. Most investigations in the 
present paper are carried out in this range. 

In the velocity field generated by DNS, we track a num- 
ber of particle pairs (1 x 10® pairs) according to the ad- 
vection equation: 



_9 
di 



(12) 



where Xi{t) is the position of the i-th particle at time 
t. We employ the 4-th order Runge-Kutta method for 
numerical integration, and the linear interpolation to ob- 
tain the velocity of each particles. Particle pairs are dis- 
tributed homogeneously with relative separation of the 
grid scale Aa; at the initial time. 



IV. RESULTS AND DISCUSSION 
A. Fixed-time statistics 

First, we briefly discuss the results obtained by stan- 
dard fixed-time statistics, which concerns the distribu- 
tion of relative separation at a certain (fixed) time. 

In Fig. 121 we plot temporal evolutions of the mean rel- 
ative separation of particle pairs, (r^), for different pow- 
ers, p = 1/2, 1, and 2. In the cases of p = 1/2 and 1, 
the generalized Richardson's relations, (r^) oc t^P^^, are 
realized, though the times at which they start and the 
time intervals in which they hold differ. These differ- 
ences result from the limited width of the inertial range 
in DNS. For example, in the case of NVll, the inertial 
range is roughly estimated as 80 < r/Ax < 320. From 
Fig. inja), it is clear that the separation PDF broadens 
rapidly, so that (r^) is contributed to by a quite broad 
range of relative separation including the dissipation, the 
inertial, and the energy containing ranges. As a result, 
the weight of each range varies with the power of the 
moment of the relative separation and thus, the range 
of time in which the contribution of the inertial range 
dominates differs accordingly. The slope of (r^) is less 
steep than Richardson's law. This is because, in the case 
oi p = 2, the contribution of particle pairs in the energy 
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TABLE II: Characteristic quantities of generated turbulent fields. Aa; represents the grid size. 



label 




Ve 


rg 


Ax 


\y 


RaA 


NVIO 


7.11 X 10"^ 


5.87 X 10"^ (0.956Aa;) 


3.44 X 10~^ 


4.50 X 10"^ 


4.54 X 10"^ 


2.50 X 10^ 


NVll 


1.44 X 10^2 


3.49 X 10"^ (1.14Ax) 


2.43 X 10-1 


3.01 X 10-^ 


3.03 X 10-2 


3.99 X 10=* 



1e+08 




t 



FIG. 2: Temporal evolution of mean relative separation 
{ritY) obtained in the case of NVll. Solid, dotted and dashed 
line refer to p = 1/2, 1, and 2, respectively. Straight lines in- 
dicate Richardson's law. The inset shows the local slope of 

(rm- 



containing range to (r^) is much larger than that in the 
cases of p = 1 and 1/2. In the energy containing range, 
relative separation process is described by the Brownian 
motion, (r^) cx t, which is less steep than Richardson's 
law. Hence, the larger the contribution of the range is, 
the less steep the slope of (r^) is. 

Even though the Richardson's law is observed for p = 
1/2 and 1, the fact doesn't necessarily support the va- 
lidity of the Richardson's law because it is reported that 
the temporal evolution of (rP) strongly depends on the 
initial separation of particle pairs |^ . Hence we have 
to check the scaling law by adopting a different method 
that is independent of the initial separation. 

In order to check self-similarity of the PDF of relative 
separation, we plot the PDF rescaled with (r^)^/^ in Fig. 
Elb). The values of {r^y/^ at t = 1.2, 1.6, 2,0, and 2.4 
are 15.3, 28.6, 47.5, and 72.0Aa;, respectively. All PDFs 
collapse well for r/(r^)^/^ < 10 and they are in good 
agreement with the similarity solution of Palm's equation 
with 6 = —0.06. However, we cannot conclude that the 
temporal evolution of the separation PDF is self-similar 
and governed by Palm's equation because Fig. ^h) is 
obtained not from separations only in the inertial range 
but in the much wider range, < r < 720Aa;. Although 
the collapse in Fig. |Sl[b) implies existence of a self-similar 
stage governed by Palm's equation, we have not had any 



reasonable explanation of the collapse. 

B. Exit-time statistics 

1. Mean exit-time 

The scale dependence of the mean exit-time obtained 
by our DNS is shown in Fig. It is observed that, al- 
though the width is narrow, there is a region consistent 
with the Bolgiano-Obukhov scaling, {Tsir; p)) oc r^/^. 
The exit-time statistics is independent of initial separa- 
tion of particle pairs if spatial scale r is large enough 
for them to forget information of their initial conditions 
[^Il0j|. In our simulation, the initial separation of particle 
pairs is Ax that is much smaller than the scales in which 
the scaling law holds. Therefore, this result indicates 
that Richardson's law is valid in the 2D-FC turbulence. 

According to the scaling law, the form of the mean 
exit-time is given by 

where C^*^'' is considered to be a universal constant. 
The inset of Fig. [S] shows a compensated plot of the 
mean exit-time. Although the weak Ra^ dependence of 
C^^'' seen in Fig. requires higher resolution of DNS to 

more accurate estimation, C^'^^ is estimated as 2.6 in 
the present DNS. 

In order to obtain statistically reliable data of exit- 
time, it is important that the residual ratio of particle 
pairs must be small enough to take in slowly separating 
particle pairs. The residual ratio at a scale r and a time 
t is defined as the ratio of particle pairs of which first 
passage times at r are less than t. If the residual ratio 
is not small enough at the time when statistics of exit- 
time are calculated, the statistics cannot reflect slowly 
separating particle pairs. Figure ^ shows the residual 
ratio of particle pairs used in the present work at the 
termination time of DNS, t = 60. It is obvious that 
the ratio is almost zero in the inertial range, so that our 
exit-time data is reliable in the sense mentioned above. 
FigureEIb) shows the scale dependence of the mean exit- 
time obtained from insufficient data and illustrates the 
importance of the residual ratio. Although the results at 
t ~ 2Q seem to have the wider inertial range than others, 
it is a fake. Therefore, to obtain statistically reliable 
data, we have to track particle pairs until the residual 
ratio becomes at least less than 0.1%. 




FIG. 3: The PDF of relative separation r obtained in the case of NVll. Different marks refer to different time: open box, 
closed box, open circle, and closed circle refer to t = 1.2, 1.6, 2.0, and 2.4, respectively, (a) the non-rescaled plot, (b) the 
rescaled plot with (r^)^^'^. The solid, dashed, and dotted lines are the similarity solutions of Palm's equation Eq. ^ with 
5 = —0.77, —1.2, and —0.06, respectively. The inset shows the linear plot of a blowup of the head region. 
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FIG. 4: Residual ratio of particle pairs. Dotted lines with 
open and closed circles refer to data obtained by DNS at reso- 
lution N = 1024 (NVIO) and = 2048 (NVll) respectively. 
The inset is the log-log plot of the ratio. 



In contrast to the inertial range, the mean exit-time 
is almost constant in the dissipation range {r/rj^ < 30 
in Fig. |3J| except in the initial transient region, because 
the relative velocity is proportional to relative separation: 
v{r) (X r. The exit-time, then, is given by 



FIG. 5: Scale dependence of the mean exit-time {Te{S;p)) 
for p — 1.1 rescaled with the dissipation time scale, tq. Open 
and closed circles refer to the results obtained by the present 
DNS at resolution N = 1024 (NVIO) and iV = 2048 (NVll), 
respectively. The inset is the compensated plot with {p" — 
l){r/ri0)^^^. The dashed line represents estimated values of 
the coefficients C^*^' ~ 2.6. 



2. PDF of exit-time and distribution of particle pairs in 
real space 



I Const. (14) 

Note that the mean exit-times rescaled by Tg collapse 
each other in Fig. [S] This is because there is only one 
characteristic time scale, Tg, in the dissipation range [2lj |. 



We plot exit-time PDF, PE{TE;r, p), in the inertial 
range in Fig. Ela), and that rescaled by the mean exit- 
time, {TE{r;p)), in Fig. Ctb). Obviously the PDFs of 
different scales collapse onto one curve in Fig[7fb). This 
means exit-time PDF is self-similar in the inertial range 
and indicates that relative dispersion process is self- 




FIG. 6: Mean exit-time and the residual ratio for several time steps (p — 1.1): open box, closed box, open triangle, and closed 
circle represent t = 20, 30, 40, and 60, respectively, (a) the residual ratio for several time step, (b) the scale dependence of the 
mean exit-time. 
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FIG. 7: The PDF of exit-time for p — 1.1 in the case of NVll. (a) the non-rescaled PDF, (b) the rescaled PDF with the mean 
exit-time {Tsir; p)). The solid, dashed, dotted and dashed-dotted lines refer to r = 80, 129, 208 and, 335 Aa;, respectively. The 
insets show a semi-log plot of the PDF. 



shuilar. 

As is clear from Fig.[7|[b), the exit-time PDF consists of 
two regions: the sharp peak and the long exponential tail. 
We call these two regions the Region-I and the Region-II, 
respectively. The qualitative difference in form between 
these indicates that the PDF reflects two different types 
of motions. In order to clarify the difference, snapshots 
of typical distribution of particle pairs in the Region-I 
and -II are shown in Figs. |S1 and El respectively, with a 
snapshot of the magnitude of T-Vorticity field \x\- 

In Fig. ISfa), on the snapshot of |x| a.i t — % {= tg), 
we plot line segments representing particle pairs of which 
exit-time T^irs', Ps) for = 80.2Aa; and Ps = 1.1 satis- 
fies 0.3 < TE/{TE{rs]Ps)) < 0.35, that is, TB(r,; p,) is in 
the shaded region in Fig. |SIb). In order to select such 



particle pairs, we first extract pairs satisfying the condi- 
tion 0.3 < Te/ {TE{rs; p)) < 0.35, and then pick out ones 
of which first passage times at are smaller than ts and 
those at prs are larger than ts- Figure Efa) is drawn by 
the same procedure. To draw Figs.|HlandO we used the 
data of NVIO. 

In 2D-FC turbulence, fine coherent structures are well 
approximated by the Burgers T- Vortex layer, so that 
shear layers are formed around the structures In 
addition, such structures are persistent in time to some 
extent. Hence, particles around the coherent structures 
are advected along them and relative separations of the 
particles expand or compress persistently. We call this 
type of motion the persistent separation. As is shown 
in Fig. |SIa), particle pairs in the Region-I appear to be 
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FIG. 8: Distribution of particle pairs in the Region-I, 0.3 < Te/{Te) < 0.35, for r = 80.2Aa; and p = 1.1 at i = 8 in the 
case of NVIO. (a) the distribution of particle pairs superimposed onto a snapshot of the modulus of T-Vorticity, \x\- Shading 
represents intensity of |x|. A particle pair is represented by a line segment (particles are located at both ends of the line 
segment). The whole computational domain is shown, (b) the PDF of exit-time rescaled with the mean. The particle pairs 
drawn in (a) belong to the gray region, 0.3 < Te/{Te) < 0.35. 




FIG. 9: Distribution of particle pairs in the Region-II, 10 < Te/{Te} < 11, for r — 80.2Ax and p = 1.1 at t = 8 in the case of 
NVIO. (a) and (b) are drawn in the same manner as Fig. |H| 
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FIG. 10: Temporal evolution of relative separation for sev- 
eral particles in the case of NVll. Different lines represent 
different particle pairs. 



along the fine coherent structures. Moreover relative sep- 
arations of particle pairs in Fig. |HIa) expand rapidly as 
can be known by their short exit-time. Therefore, Fig.|Sl 
supports the picture of persistent separations. 

In contrast to the Region-I, the distribution of particle 
pairs in the Region-II is in disorder (see Fig.EIa)). Be- 
cause the particle pairs contained in the Region-II have 
long exit-time, the longitudinal relative velocity of them 
remains small positive during the passage from to pvs 
or moving direction of relative separation of the pairs 
changes from expansion to compression at least once be- 
fore they reach pr^. However it is quite unlikely that 
relative velocity remains small positive for a long time, 
so that the latter is probably the main mechanism of 
the formation of the Region-II. Hence, it is expected 
that both of expanding and compressing pairs are con- 
tained in Fig.l^Ja). In fact, there are several pairs placed 
along coherent structures, which probably are expand- 
ing. Note that their relative separations are not confined 
to be larger than r^. In addition, there are also pairs 
across the structures, which are presumably compressing 
because the structures are approximated by the Burg- 
ers T- Vortex layer. Moreover, some pairs are located at 
positions where the fine coherent structures twist. The 
twisted regions are generated when well-stretched struc- 
tures lose their activities and are folded, or when the 
structures are generated by plumes. Therefore the par- 
ticle pairs in such regions are changing their moving di- 
rections from expansion to compression or the opposite. 

Figure EH shows five typical evolutions of particle-pair 
separation. Sample 3 expands persistently without any 
strong compression but Sample 4 experiences both per- 
sistent expansion and compression several times. The 
former corresponds to motions along structures and be- 
longs to the Region-I; the latter does across structures or 
twisted region and belongs to the Region-II. Note that 



FIG. 11: Estimated values of 5 in the case of NVll. Dotted 
lines with open box, closed box, open circle, closed circle, open 
triangle, closed triangle, and reverse open triangle denote p = 
l.l^/^ 1.1^/*, l.l^/^ 1.1, l.l^ l.l", and 1.1*, respectively. 
Solid and dashed lines indicate the mean value, 5 = —0.77, 
and the value of the Richardson's case, & — —1.2, respectively. 
The shaded region represents the standard deviation. 



even a single evolution contains motions that are catego- 
rized into the Region-I or into the Region-II. These facts 
are consistent with the assumptions of the self-similar 
telegraph model: relative separation process consists of 
persistent expansion and compression with some random 
transition mechanisms. 



3. Estimation of 5 

It is expected that expansion and compression of rela- 
tive separation differ in persistency; the former is experi- 
enced mainly by particle pairs along coherent structures 
and the latter does by those across the structures. Since 
the auto-correlations of strain and temperature along the 
structures have longer characteristic lengths than those 
across them in 2D-FC turbulence expanding mo- 

tions must be more persistent than compressing ones. So 
as to confirm this consideration, we investigate the slope 
of the tail of the exit-time PDF. Then, with Eq. (0), we 
estimate 5 that describes the difference in persistency be- 
tween expansion and compression of relative separation. 

Note that using the asymptotic form of the exit-time 
PDF calculated from Palm's equation, Eq. (0), may be 
justified as follows. Introducing a time scale t, the l.h.s. 
of the self-similar telegraph model, Eq. Q , is rewritten 
as follows: 



l.h.s. 



Tc{r) 1 d'^P 1 dP 



(15) 



where t = t/t. In the case that t ^ Tc{r)/\ at a certain 
spatial scale r, the first term of the l.h.s. can be ne- 
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glected, and thus, the model is reduced to Palm's equa- 
tion. Hence, the tail of the exit-time PDF by the tele- 
graph model, which consists of slowly-separating particle 
pairs, must coincide with that by Palm's equation. It is 
also easy to show that the mean exit-time is the same 
between Eqs. (0) and ©• We, therefore, adopt Eq. © 
for the estimation of S. 

Figure^]shows the estimated values of 6 in the incrtial 
range. The mean value of 5 is —0.77 and the standard 
deviation is ±0.17. This negative value denotes that ex- 
panding motions are more persistent than compressing 
ones, which supports the above expectation. 

In addition, S is larger than that of the Richardson case 
(the zero-drift), S — —1.2 = Sq. Because (5 = A+ — A~, 
S > Sq means A~ < X'^ — Sg = Xq , that is, the com- 
pressing motion of relative separation is more persistent 
than that of the Richardson case. This fact results in the 
negative drift in the self-similar telegraph model, Eq. ||2J). 
The negativity of the drift can be accepted considering 
that coherent structures in 2D-FC turbulence are string- 
like and scale-transversal ones That is, because the 
coherent structures in 2D-FC turbulence are not nested 
in scales, there are few obstacles to compressing motions. 
On the other hand, the typical structures in 2D-IC tur- 
bulence are nested vortices called "cat's eye in a cat's 
eye". A separation process of particle pairs in 2D-IC 
turbulence, thus, is a step-by-step one each step of which 
consists of a trapping by one of nested vortices and sud- 
den separation into a next larger vortex |lOj | ; compressing 
motions are probably blocked by the nested vortices. In 
fact, the estimated values of 6 from DNS of 2D-IC tur- 
bulence by Boffetta & Sokolov 's'l and Goto & Vassilicos 
[Toj are smaller than that of the Richardson case J_6J ; the 
compressing motion of relative separation is less persis- 
tent than the Richardson case. This results in the pos- 
itive drift in Eq. (0). We, therefore, conclude that the 
drift term in the self-similar telegraph model, Eq. ||2Jl, 
reflects the characteristics of coherent structures of the 
flow. 



4- Comparison of the Region-II with a solution of Palm's 
equation 

In Fig. El we compare the PDF of exit-time obtained 
by DNS with that calculated from Palm's equation by 
Eq. ||B4|). In the Region-I, the form of the PDF by DNS 
is totally different from that by Palm's equation. In the 
Region-II, however, the two PDFs collapse onto a single 
curve even when p is very small, and thus, the separa- 
tion process of particle pairs in the Region-II can be de- 
scribed by Palm's equation, Eq. ©. This indicates that 
the motion of particle pair is diffusive in the Region-II; 
the separating process from rp~^ to r doesn't affect that 
from r to rp. It is also observed that the form of the 
exit-time PDF varies depending on the value of p in Fig. 

The larger the value of p is, the larger the proportion 
of the Region-II becomes. If p is sufficiently large, the 
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FIG. 12: Rescaled PDF of exit-time (r = lOTAx) in the case 
of NVll and the prediction of Palm's equation {S = —0.77) for 
several values of p. Solid and dashed lines refer to the results 
of DNS and that obtained from Palm's equation through the 
relation Eq. <B4ll . 

whole PDF seems to be occupied by the Region-II and 
the relative separation process is substantially described 
by Palm's equation, Eq. ©. 

As mentioned in the previous subsection, the tail of the 
exit-time PDF calculated from the self-similar telegraph 
model, PgT"^ (Te ; r, p) , must agree with that calculated 
from Palm's equation, P^\TE',r, p). In addition, as the 
value of p gets larger, the region of pj^^TE^r, p) over- 
lapping with p'^\Te\ t, p) expands. This is because, for 
large values of p, most of exit-times of particle pairs are 
longer than the characteristic time Tc{r)/X. Hence, the 
self-similar telegraph model probably has the same prop- 
erties as the results shown in Fig. 1121 

5. Characteristics of the Region-I and estimation of A 

Figure [T^ shows the exit-time PDF for NVll rescaled 
with {TE{r; p)) and normalized by their peak values. It 
is clear that the shape of the Region-I is independent of 
p. That is, there exists a p-indepcndcnt function repre- 
senting the shape of the exit-time PDF in the Rcgion-I, 
Pe\Te)- Here, ~ denotes rescaling by the mean exit- 
time, (Tsir; p)). We assume that Pe is normalized as 
PE{fE)dfE = 1. Then, in the Region-I, the exit- 
time PDF is written as 

PE{fE;p)^W{p)P^j^\fE), (16) 

where W{p) is a normalization factor depending on p. 
The exit-time PDF, Pe{Te;p), takes the peak value at 
ff'^'^. It is estimated that Tf'^'^ w 0.3 in Fig.[T3| which 
can be regarded as a characteristic time scale of the 
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telegraph model is given by 

^ A 1 



As{2s~5y 
Then, A is calculated as follows: 

2s -5 

A - 



^(i) 



(20) 



(21) 



If we assume that T, 



(I) 



E 



T^max 



^ 0.3, then, A is estimated 
as 5.2. However, the value of T^'^ is not necessarily well 
defined. Roughly speaking, it may take a value satisfying 
0.19 < r|?^ < 0.45, if assuming f|?^ is in the half-value 
width of the Region-I around T]^^^. Then, the estimated 
value of A is in the range, 3.5 ;< A < 8.3. 



FIG. 13: Rescaled PDF of exit-time (r = lOTAx) in the case 
of NVll. The peak value is normalized to be unity. Different 
lines represent different values of p: p = 4.59, 2.14, 1.46, 1.21, 
1.1, 1.05, 1.02, and 1.01 from the curve with the lowest cut-off 
rescaled exit-time to that with the highest. The inset is the 
same plot in linear scale. 



Region- I. As shown in the inset of Fig. US P^e\Te) dis- 
tributes sharply around T^^^^. 

If the Rcgion-I is formed by particle pairs of which rel- 
ative separations expand persistently according to v{r) 

by the pass through time from r to pr 



the exit-time in the Region-I, T^\r;p), is given 



pr 



dr 
v{r) 



(p" - ly 
71 



(17) 



Then, using the general form of the mean exit-time, Eq. 
(|B2ll . the exit-time rescaled by {TE{r; p)) in the Region-I, 



Ti\r;p), is 



Ti\r;p) 



Ti\r;p) 1 
{TE{r-p)) sCeA' 



(18) 



which is independent of p. Thus, P^\Te) can be con- 
nected to the PDF of A, P^{A), which is the PDF of 
Lagrangian relative velocity, w(r), rescaled with r^~'^: 



sCeA 



= sCeA^P^{A). (19) 



Note that P^iA) is expected to be independent both of 
r and p. Since, in the self-similar telegraph model, the 
distribution of the coefficient. A, of the relative velocity 
is not considered, Pa{A) is a ^-function: Pa{A) = 5{A — 
Ac). Accordingly, the Region-I of the exit-time PDF is 
approximated by a ^-function at T^^' in the model. 

Combining Eqs. © and (|18|l . we can estimate A from 
f^^^. From Eq. ®, Cij in the case of the self-similar 



V. CONCLUDING REMARKS 

We have investigated relative dispersion in 2D free con- 
vection turbulence by direct numerical simulation. In 
the inertial range, where the entropy cascade dominates, 
we have confirmed with exit-time statistics that relative 
dispersion satisfies the Bolgiano-Obukhov scaling and, 
therefore, is self-similar. It was also shown that the 
exit-time PDF, PE{TE',r^p), is divided into two parts, 
the Region-I and -II, and that both of them satisfy the 
scaling-law and the self-similarity. Pe{Te', r, p) is written 
as the following form: 

Pe{Te; r, p) = H{Ti\r- p) - TE)P'i^ [TeW, p) 

+ HiTE-Tt{r;p))Pi'\TE;r,p)., (22) 

where Tf^(r;p) is a division time-scale (exit-time) be- 
tween the Region-I and -II, and H{x) is a smoothed step 
function such that H{x) = if a; ^ and H{x) — 1 if 

X 3> 0. pj^\TE;r, p) and P^^^\TE\r, p) correspond to 
the exit-time PDF of the Region-I and -II, respectively. 
The investigation of the distribution of particle pairs in 
the real space indicates that the Region-I and -II are 
formed, respectively, by particle pairs expanding along 
coherent structures and by those experiencing turns be- 
tween expansion and compression (Figs. |S1 and l^]). 

Figures El and El show the following characteris- 
tics of P^^'(TB;r,p) and P];"^ (T^ ; r, p) . The form of 

Pe^ {Te ; r, p) is independent of p if it is rescaled with 
the mean exit-time. Moreover, if we assume that the 
Region-I is constituted by persistently-separating parti- 
cle pairs, the form of p],^^ (Te; r, p) is related to the PDF 
of Lagrangian relative velocity, v{r) — Ar^"". That is. 



P'i\TE;p) 



Pa 



1 



W{p) 

sCeTI " \sCeTe 



(23) 



where W{p) is a p-dependent normalization factor, and 
P^{A) is the PDF of v{r) rescaled with r^~*. On the 
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the essential characteristics of the relative separation pro- 
cess, the distribution of the relative velocity is crucial to 
understand further details of separation processes and 
their relation to the coherent structures, as suggested by 
the exit-time PDF obtained by DNS. In fact, we have ob- 
tained preliminary results showing that the distribution 
of the relative velocity in the inertial range is tightly con- 
nected to that in the dissipation range ^] . This strongly 
indicates that the distribution directly relates to the co- 
herent structures because their length are of the order of 
the inertial range and their width of the order of the dis- 
sipation range |l8j | . An extension of the model to include 
the distribution of relative velocity, as well as an exten- 
sion to the dissipation range, will be done in the future 
work. 



FIG. 14: Similarity solution (fixed-time PDF of relative sepa- 
ration) of the self-similar telegraph model in the case of 2D-FC 
turbulence. The parameters used in this figure are estimated 
from DNS data for NVU by using the exit-time PDF. Dif- 
ferent lines refer to different values of A. \ = oo corresponds 
to the similarity solution of Palm's equation. The inset is a 
blowup of the head region in the linear scale. 

other hand, P^^\Te', r, p) agrees with the exit-time PDF 

calculated from Palm's equation, p'^\Te] r, p), the form 
of which varies depending on p even when it is rescaled 
with {TE{r;p)). Hence Pe{Te] p) is written as 

PE{fE-.p) = H{fi^{p)~fE)^^^PA f^^) 

sCeTI \sCeTeJ 
+ H{fE~ft{p))PP{fE;p). (24) 

These results support the self-similar telegraph model. 
In the model, Pe\Te) is approximated by a (5-function. 
On the other hand, for slowly-separating particle pairs, 
the model is approximated by Palm's equation. That is, 
the head region of the fixed-time PDF as well as the tail of 
the exit-time PDF is approximately described by Palm's 
equation. By taking advantage of these characteristics of 
the exit-time PDF, we can estimate the parameters of the 
model, A and d' (S): X is estimated in the Region-I and S is 
in the Region-II. By our DNS of 2D-FC turbulence, these 
parameters are estimated as 3.5 ;< A < 8.3 and ct « 0.43 
{S « —0.77). The estimated value of a corresponds to 
the negative drift in the model. This is the opposite to 
the 2D-IC turbulence case T^. We speculate that the 
difference in the sign of a is caused by the difference 
in coherent structures. With these parameters, we can 
numerically calculate the similarity solution of the self- 
similar telegraph model, which is shown in Fig. 1141 1^. 
However we cannot compare the solution directly with 
the results of DNS because the inertial range achieved 
by the DNS is too narrow. We need DNS with higher 
resolutions to the comparison. 

Although the self-similar telegraph model can capture 
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APPENDIX A: PROPERTIES OF 2D-FC 
TURBULENCE 

The governing equations of the 2D-FC turbulence are 

V-M = 0, (Al) 

du , , Vp . , . , 

— (m • V)m = + vl\u - agTeg, (A2) 

ot po 

^ + {u-V)T = kAT, (A3) 

where u, T, and p represent velocity, temperature, and 
pressure field, respectively. Eg is the unit vector in the 
direction of the gravity, v, k, po, a, and g are the kine- 
matic viscosity, the molecular diffusivity, the mean den- 
sity of the fluid, the thermal expansion coefficient, and 
the gravitational acceleration, respectively. 

The 2D-FC turbulence has two important proper- 
ties: the Bolgiano-Obukhov scaling and the fine coherent 
structures. 

1. Bolgiano-Obukhov scaling 

In this system, the integral of the squared temperature, 

S=JJ,j^d-. (A4) 

is a conserved quantity in the ideal case, where is the 
volume of the system. We call this quantity entropy for 
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convenience. The entropy cascades from large to small 
scales similar to the energy cascade in the 3D-NS turbu- 
lence, which leads scaling laws of the energy and entropy 
spectra, E{k) and S{k): 

E{k) - KEee^'\agf"'k-^^'\ (A5) 
S{k)^Ksee^"'{ag)-^/'>k-''\ (A6) 

where Ke and Ks are considered to be universal con- 
stants, and te is the entropy dissipation rate defined as 
follows: 



VT ■ VTdxj 



(A7) 



^-1/4^5/8 



In this paper, (•) denotes an ensemble average. These 
scaling laws, Eqs. lASp and ljA6p . are called the Bolgiano- 
Obukhov scaling _2j • 

Entropy dissipation length scale rjg and time scale Tg 
is estimated by dimensional analysis with eg, ag, and k; 

(A8) 
(A9) 

These scales correspond to the Kolmogorov scales in the 
3D-NS turbulence. There are several other quantities 
characterizing 2D-FC turbulence: the thermal Taylor mi- 
croscale A, and the Rayleigh number at A, Ra^. These 



Te = ee ^^^{(^gY 



are defined as A 
where 



^(Ax + Aj,) andRaA = i(RaA, -f RaA J, 



Aa; 



A„ = 



{(dT/dx) 



RaA, 



RaA„ 



/ {{dT/dyY)' 

KV ' 

ag{TY'^\l 



(AlO) 

(All) 

(A12) 
(A13) 



According to the Bolgiano-Obukhov scaling, Richard- 
son's law in the 3D-NS turbulence is modified as follows: 



(r(if)=g,e//2(«g)Pi5p/2 



(A14) 



where r{t) is relative separation of passive particles at 
time t and gp is considered to be a universal constant. In 
addition, Richardson's diffusion equation of the separa- 
tion PDF P{r,t) is modified as follows: 



dP 
'dt 



d_ 

dr 



{koe^iagr^'r^^^jr 



dr 



(A15) 



where P{r, t) is the probability density of relative sepa- 
rations at r and t, and fco is considered to be a universal 
constant. This equation has a self similar solution if the 
initial condition is the delta function. For P(r, 0) — S{r), 



P{r,t) 



2TrrCB. 



ee{agY{koty 



■ exp 
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r.2/5 



4 koey'{ag)VH^ 



(M6) 



where Cr is the normalization constant. 



2. Fine coherent structures 

We introduce a vector quantity called T-Vorticity, 
x(a;, t), as follows: 



dT dT\ 
dy' dx J 



(A17) 



Figure E| shows a snapshot of the magnitude of the 
T-Vorticity and the strain field. It is observed that there 
are linearly concentrated area of T-vorticity. Its width 
and length are an order of the entropy dissipation scale 
and the integral scale, respectively [l^- We call these 
fine coherent structures. At the same location of the 
structures, we can observe linearly straining regions in 
Fig. cab). 

T-Vorticity has similar properties to vorticity in the 
3D-NS turbulence. The evolution equation of the T- 
Vorticity has the stretching term and is the same as that 
of vorticity in 3D-NS system: 



(A18) 



Accordingly, there is a solution corresponding to the 
Burgers vortex in the 3D-NS system, which is called the 
Burgers T-vortex layer. It is numerically confirmed that 
fine coherent structures in the 2D-FC turbulence are well 
approximated by the Burgers T-vortex layer [l^ . This is 
similar to fine coherent structures such as worms in the 
3D-NS turbulence which are well approximated by the 
Burgers vortex. 



APPENDIX B: EXIT-TIME STATISTICS 

The exit-time statistics is one of the scale-fixed statis- 
tics. The exit-time, Te{R] p), is defined as 



TE{R;p)=TFipR)-TFiR), 



(Bl) 



where Tf(R) is the time when a relative separation r{t) 
reach a threshold R for the first time (first-passage time) 
Hi ^1 HSl ■ According to the scaling law of the character- 
istic time, Tc cx r®, and the additivity of the mean, the 
form of the mean exit-time is expected as follows: 



{TEiR;p))=CE{p'-iy, 



(B2) 



where Ce is a constant depending on the system. 
In the case of the Bolgiano-Obukhov scaling, Ce — 



-1/5 



where C 



(BO) 



is considered to be a 



non-dimensional universal constant. 

The exit-time statistics has two advantages over usual 
fixed-time statistics: 

1. The exit-time statistics can specify a spatial scale 
by choosing a threshold R. Therefore we can ex- 
tract information of the inertial range if both R and 
pR are in the inertial range. 




FIG. 15: A snapshot of T-Vorticity and strain field, (a) a snapshot of T-Vorticity field, (b) that of strain field at the same 
location and time. Shading represents intensity of the field. 



2. The interval between thresholds (the width for av- 
eraging) can be controlled by p. This means that 
we can control the degree of coarse graining of dis- 
persion process. 

In order to calculate exit-times, we prepare a set of 
thresholds i?„ = p^Ro, where n is a positive integer, and 
record the time at which r(t) reach Rn for the first time 
for every particle pairs and every n. In the present work, 
we set p to 1.1^/^ (w 1.01) and Rq to the grid size Ax. 

If Palm's equation, Eq. ©, describes relative disper- 
sion well, from the equation we can calculate the proba- 
bility density function (PDF) of exit-time PeITe; R, p), 
which represents the probability density for a particle 
pair of which exit-time from R to pR is Tg, Q. In the 
present case, to calculate the PDF, first we solve Eq. Q 



with the initial condition 

P(r,0) ^d{r-R), 



(B3) 



where the boundary conditions are the reflecting condi- 
tion at r = and the absorbing condition at r = pR. 
And then, the PDF is obtained as the time derivative of 
the probability of r{t) < pR: 



d f 

Pe{Te;R,p) ^--n P{ 

"'^ Jr<pR 



r, t)dr 



(B4) 



t=TE 



We have numerically calculated the PDF of exit-time by 
the above procedure and compared with the PDF ob- 
tained by our DNS. 
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